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REVIEW OF EXPERIMENTAL INVESTIGATIONS OF LIQUID-METAL HEAT TRANSFER 1 


By Bernard Lubarsky and Samuel J. Kaufman 


SUMMARY 

The experimentally obtained results of various investigators 
of liquid-metal heat-transjer characteristics were examined and 
fund to be not always directly comparable because of differences 
in experimental apparatus or in methods of calculation. The 
experimental data were therefore reevaluated using assumptions 
and methods as consistent as possible and then compared with 
each other and with theoretical results. 

The reevaluated data for both local fully developed and average 
Nusselt numbers in the turbulent-flow region were still found to 
have considerable spread, with the bulk of the data being lower 
than predicted by existing anaylses. An equation based on 
empirical grounds which best represents most of the fully 
developed heat-transfer data is 

Nu= 0.625 Pe 0i 

where Nu represents the Nusselt number and Pe, the Peclet 
number. The theoretical prediction of the heat transfer in the 
entrance region was found to give lower values, in most cases, 
than those found, in the experimental work. 

The theoretical and experimental results for the ratio of local 
Nusselt number to fully developed Nusselt number were inte- 
grated to obtain predictions for the ratio of average Nusselt 
number to fully developed Nusselt number for a range of Peclet 
numbers and length-diameter ratios. Most of the experimental 
data fall between 60 to 80 percent of the predicted values. 

The experimental evidence was insufficient to serve as a basis 
for any conclusion concerning liquid-metal heat transfer in the 
laminar or transition flow regions. 

INTRODUCTION 

The use of liquid metals as heat-transfer media is presently 
of considerable interest. A number of theoretical and experi- 
mental investigations to determine the heat-transfer charac- 
teristics of liquid metals have been made by various 
investigators (refs. 1 to 26). In the literature, the results of 
the experimental investigations often have been compared 
with each other and with the results of theoretical investiga- 
tions. During the course of investigations of liquid-metal 
heat-transfer characteristics at the NACA Lewis laboratory, 
the work of the various experimental investigators was 
carefully examined. It was found that different investiga- 
tions were not always directly comparable because of differ- 


ences in the experimental apparatus or in the methods of 
calculation. Some of the differences found were: 

(1) Liquid-metal physical properties that differed from 
those currently accepted were sometimes used. 

(2) At times, centerline temperatures in and out of the 
test section were measured rather than “mixing-cup” 
temperatures. 

(3) Some of the experiments were conducted with uni- 
form heat input to the wall of the test section, while others 
more closely approached constant wall temperature. 

(4) Some investigators measured the combined heat- 
transfer coefficient in a tube and concentric annulus; different 
methods were used to obtain the individual coefficients. 

(5) Some investigators measured local fully developed 
heat-transfer coefficients; others measured average over-all 
coefficients. 

(6) The velocity profiles entering the test section varied; 
some approached a fully developed turbulent profile, while 
others were more nearly uniform. 

(7) Different length-diameter ratios of the test section 
were used. 

The differences in experimental apparatus of items (6) 
and (7) affect only the average heat-transfer coefficient and 
not the fully developed coefficient. 

Because of the differences in experimental apparatus and 
methods of calculation listed, the experimental data of 
references 1 to 26 were reevaluated using consistent assump- 
tions and methods in order to permit a better intercom- 
pax ison of the experimental results and comparison with the 
results of theoretical investigations. 

SYMBOLS 

a constant 

c p specific heat, Btu/(lb)(°F) 

D equivalent or hydraulic diameter, ft 

Di annulus inner diameter, ft 

D 0 annulus outer diameter, ft 

f friction factor 

G weight flow per unit area, lb/(hr)(sq ft) 

Gz Graetz number, PeD/l or PeT)jx 

k thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 

l length of test section, ft 

m constant, eq. (10) 

Nu Nusselt number, UD/k 

n constant, eq. (10) 


1 Supersedes NACA TN 3336, “Review of Experimental Investigations of Liquid-Metal Heat Transfer,” by Bernard Lubarsky and Samuel J. Kaufman, 1955. 
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Pe Peclet number, RePr, GDc P /k 

Pr Prandtl number, c r fi/k 

Re Reynolds number, GD/ti 

St Stanton number, U/c p G 

t c fluid centerline temperature, °F 

f„, fluid bulk temperature, °F (“bulk temperature” as 
used in this report is synonymous with “mixing-cup 
temperature” and “mixed mean temperature”) 
t„ wall temperature, °F 

U heat-transfer coefficient, Btu/(lir)(sq ft)(°F) 

x distance along test section, ft 

jj, fluid bulk viscosity, lb/ (hr) (ft) 

Subscripts: 
an annulus 

av average 

/ fully developed 

x at station x 

PROCEDURE 

The experimental data of the various references were 
reevaluated as consistently as possible, plotted as Nusselt 
number against Peclet or Graetz number or against both, and 
the results compared with theoretical predictions. These 
three steps will be discussed in reverse order, because some 
of the methods used in reevaluating the data were deter- 
mined by theoretical considerations. 

THEORETICAL INVESTIGATIONS OF LIQUID-METAL HEAT TRANSFER 

The following discussion gives a brief description of some of 
the results of theoretical investigations and is not intended 



Peclet number, Pe 


Figure 1. — Theoretical predictions of fully developed Nusselt numbers 
for heat transfer to liquid metals in turbulent flow in round tubes. 


to be complete. All the theoretical investigations discussed 
consider only the turbulent-flow region. 

Fully developed heat-transfer coefficients. — Heat-transfer 
coefficients for liquid metals in turbulent flow with fully 
developed velocity and temperature profiles have been pre- 
dicted by a number of investigators using somewhat different 
assumptions: 

(1) Uniform heat input to the wall; round tubes: The most 
frequently analyzed case is that of heat transfer to a round 
tube with uniform rate of heat input along the length of the 
tube. This case was investigated by Martinelli (ref. 27) 
using the “momentum transfer analogy.” Lyon (ref. 6) 
found a simplified equation which approximated Martinelli’s 
more complex relation. This equation, which is recom- 
mended by the Liquid-Metals Handbook (ref. 28), is 

Nu f =7.0-\-0.025 Pe? 8 (1) 

Cope (ref. 29) investigated the possibility of assuming that 
the modified vorticity transfer analogy applied to the tur- 
bulent core of the fluid, while the momentum transfer 
analogy applied to the boundary layer and buffer layer. 
Kennison (ref. 30) assumed that the heat transfer is analogous 
to the transfer of vorticity for turbulent fluid flow in a long 
straight pipe. Deissler (ref. 31) modified the momentum 
transfer analogy to allow for heat transferred by conduction 
to or from a turbulent particle as it moves radially in the 
tube. Deissler’s analysis is for a Prandtl number of 0.01. 

Some of the results of these various investigations are 
shown in figure 1. The experimental results for fully 
developed heat transfer in a round tube with uniform heat 
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Figure 2. — Theoretical predictions of fully developed Nusselt num- 
bers for heat transfer to liquid metals in turbulent flow in annuli 
and between flat plates. 
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input will be compared with Lyon’s equation (eq. (1)) inas- 
much as this is the equation recommended by the Liquid- 
Metals Handbook and most commonly used in practice. 

(2) Uniform wall temperature; round tubes: The fully 
developed heat-transfer coefficient in turbulent flow in a 
round tube with a uniform wall temperature has been 
investigated by Seban and Sliimazaki (ref. 32) using the 
momentum transfer analogy; they give, as an approximate 
relation, the equation 

iVw / =5.0+0.025 Pe r 0 - 8 (2) 

This equation is also plotted in figure 1. The Liquid-Metals 
Handbook lists the equation as 

iV%=4.8 + 0.025 Pe t o a (3) 

and gives the work of Seban and Sliimazaki as a reference. 
The experimental results for fully developed heat transfer 
in a round tube with a uniform wall temperature will be com- 
pared with Seban and Shimazaki’s equation (eq. (2)). 

(3) Uniform heat input; annuli: Very little theoretical 
work has been done on the fully developed heat-transfer 
coefficient in annuli. For thin annuli (diameter ratio <1.4) 
the Liquid-Metals Handbook recommends the use of the 
theoretical relation proposed by Seban (ref. 33) for heat 
transfer to parallel plates with heat through one side only: 

A%/,„k=5.8+0.020 Pe r ,an 8 for /+/£>,< 1.4 '"Jg (4) 

For annuli of diameter ratio greater than 1.4, the Liquid- 
Metals Handbook lists an equation which approximates the 
results of Bailey (ref. 34) and is of the form suggested by 
Werner, King, and Tidball (ref. 7): 

WM / , a „=0.75(H o /H i ) 0 - 3 (7.0 + 0.025 Pq„ 0 - 8 ) iovDJD^lA 

(5) 

Equations (4) and (5) are plotted in figure 2. The experi- 
mental data on heat transfer in annuli will be compared with 
these equations. 

Local heat-transfer coefficients in entrance region. — Heat- 
transfer coefficients in the entrance region have been calcu- 
lated by several investigators for a number of different cases. 
Poppendiek, Palmer, and Harrison (refs. 26, 35, and 36) have 
analyzed the case of uniform wall temperature for various 
different entering velocity profiles; the analysis assumes that 
the eddy diffusivity of heat is negligible when compared with 
the molecular diffusivity and consequently is intended only 
for low Reynolds numbers. The analysis is independent of 
Prandtl number. Deissler (ref. 37) analyzed the case of 
uniform heat input at the wall, with a fully developed velocity 
profile at the entrance; the numerical calculations were 
carried out only for a Prandtl number of 0.01. Seban and 
Shimazaki (ref. 38) have made calculations for the case of 
uniform wall temperature and fully developed velocity profile 
at the entrance for a Prandtl number of 0.01 and Reynolds 



(a) Ratio of entrance Nusselt number to fully developed Nusselt 
number against Graetz number. 
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(b) Nusselt number against Peclet number. 

Figure 3. — Theoretical predictions of heat transfer to liquid metals 
in turbulent flow in round tubes. 

numbers of 10 4 and 10 5 . The results of the analyses of 
Poppendiek and Palmer and of Deissler are shown in figure 3. 

Average heat-transfer coefficients. — Predictions of average 
heat-transfer coefficients can be made by integrating the 
predictions for local heat-transfer coefficients over the length- 
diameter ratio of the tube in question. Heat transfer in the 
entrance region, however, has been analyzed for only rela- 
tively specialized cases. Therefore, the experimental results 
for average heat-transfer coefficients will first be compared 
with equations (1) and (2), even though these equations are 
derived for fully developed heat-transfer coefficients. Later 
in the report, a comparison will be made with the average 
heat-transfer coefficients on the basis of the analytical 
evidence. 
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Figure 4. — Theoretical predictions of Martinelli (ref. 27) for fully 
developed temperature profiles for heat transfer to liquid metals in 
round tubes. Prandtl number, 0.022. 


Temperature distribution. — The fully developed temper- 
ature distribution due to heat transfer to a liquid metal in 
turbulent flow in a round tube has been predicted on theoret- 
ical grounds by several investigators. The predictions of 
Martinelli (ref. 27) are shown in figure 4 for a Prandtl num- 
ber of 0.022. Martinelli, using his own predicted values for 
the temperature distribution, calculated the ratio of the 
temperature differences (t w —t m )/(t w —t c ) as a function of 
Reynolds and Prandtl numbers. Martinelli’s results are 
shown in figure 5(a) for Prandtl numbers pertinent to liquid 
metals. Martinelli also calculated values of ( t w — t m )/ 



(a) Mercury. 

Figure 6. — Physical properties of liquid metals (ref. 28). 


{t w — t-c) for fully developed flow between flat plates with heat 
flow through both walls with uniform heat flux. These 
results are shown in figure 5(b). 






Figure 5. — Theoretical predictions of Martinelli (ref. 27) of ratio — for heat transfer to liquid metals. 
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The heat-transfer parameters were evaluated using the 
same method of calculation for each individual reference as 
was used by the authors of that particular reference, with the 
following exceptions: 

(1) All physical properties of liquid metals were taken 
from the second edition of the Liquid-Metals Handbook (ref. 
28). These properties are shown in figure 6. 

(2) When an investigator measured the combined liquid- 
metal heat-transfer coefficient in a tube and concentric 
annulus, the individual heat-transfer coefficients were ob- 
tained by assuming that the ratio of the Nusselt number in 
the tube to the Nusselt number in the annulus is determined 
by equations (1), (4), and (5): 


Nu 

Nu an 


Nu 7.0+0.025P<? 0 - 8 
Nu an 5.8+0.020-Pea,, 0 - 8 

7.0+0.025Pe 0 - 8 

= 0.75(H„/H f )°- 3 (7.0+0.025Pe 

.01 2 


for #^<1.4 (6) 

. for D 0 /D t >1.4 (7) 


.008 


— 01 
"O -O L 

g § £ 


.004 


O QJ 



m .30 



Temperature, °F 

(c) Lead-bismuth eutectic. (d) Sodium. 

Figure 6. — Concluded. Physical properties of liquid metals (ref. 28). 
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Because of the lack of theoretical work on average heat- 
transfer coefficients, particularly in annuli, the same ratios 
which have been assumed for the fully developed Nusselt 
numbers will be assumed for the average Nusselt numbers. 

It is important to note that in most of those tests in which 
the combined coefficient in a tube and concentric annulus 
was measured, the Reynolds number in the annulus was 
smaller than the Reynolds number in the tube. Quite often 
the flow in the annulus was in the transition flow region, 
while the flow in the tube was in the turbulent-flow region. 
Inasmuch as there are no predictions for liquid-metal heat 
transfer in the transition region, equations (6) and (7) will 
be used to separate the tube and annulus heat-transfer 
coefficients even when the flow in the annulus is in the 
transition region. This procedure is open to question, and 
the interpretation of the data calculated by this procedure 
may be inaccurate. 

(3) In those tests in which the centerline temperature of 
the fluid was measured instead of the bulk temperature, the 
temperature difference between the wall and the bulk fluid 
will be calculated from Martinelli’s relation for ( t w —t m )/ 
(; t w —t c ) (fig. 5). Martinelli’s prediction of 
for flat plates with heat flowing through both sides will be 
used for annuli inasmuch as no other predictions covering as 
broad range of Reynolds and Prandtl numbers are available. 

REEVALUATION OF EXPERIMENTAL DATA 

The experimental investigations of references 1 to 26 will 
first be discussed individually and then compared with each 
other and with theoretical investigations. 

The experimental work of the various investigators will be 
discussed in a chronological order determined by the publica- 
tion date of the original manuscript. 

Styrikovich. and Semenovker. — Styrikovich and Semenov- 
ker (ref. 1) investigated heat transfer to mercury as part of 
their investigation of the mercury-steam binary power cycle. 
They used a series of five tubes for test sections, each about 
106 inches in length with 0.63-, 0.87-, 1.58-, 1.67-, and 1.97- 
inch diameters. The tubes were heated by external electric 
heaters. Thermocouples were placed 17.2 inches apart on 
the outside surface of each tube. The bulk fluid tempera- 
ture in the test section was calculated by adding to the inlet 
temperature the temperature rise corresponding to the heat 



Figure 7. — Comparison of variations of Prandtl number of mercury 
with temperature of Liquid-VIetals Handbook (ref. 28) and of 
Styrikovich and Semenovker (ref. 1). 
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Figure 8. — Reevaluated data of Styrikovich and Semenovker (ref. 1) 
for fully developed heat transfer to mercury in round tubes. 

input. The velocity profile of the mercury entering the test 
section was essentially fully developed. The method of 
heating the mercury approximated uniform heat input to the 
wall. The heat-transfer coefficient was calculated for only 
the central portion of the tube. The coefficients presented 
are essentially the fully developed heat-transfer coefficients. 

The physical properties used in evaluating the heat- 
transfer coefficients are not listed, but the Prandtl number is 
tabulated over a range of temperature from 32° to 1112° F. 

These Prandtl numbers are lower than the values in reference 
28, which lists values of Prandtl number for temperatures 
up to 600° F. The values of Prandtl number of Styrikovich 
and Semenovker and of reference 28 are shown in figure 7. 

Since the specific heat and viscosity in the temperature range 
used are essentially the same in reference 28 as those reported 
in the International Critical Tables (1929 edition), the 
inaccuracies in Prandtl number may be assumed due to 
incorrect values of thermal conductivity. It appears that 
Styrikovich and Semenovker used the thermal-conductivity 
data of Gelhoff and Neumeier, which have been found to be 
high (ref. 12). It was deemed advisable to recalculate the 
data of Stjn-ikovich and Semenovker using the values of 
thermal conductivity from reference 28. The precise tem- 
perature level of the various data points is not reported, 
but the average temperature level is given as about 932° F. 

At this temperature, Styrikovich and Semenovker list a 
Prandtl number of 0.0056. Reference 28 presents Prandtl 
number data up to 600° F which when extrapolated to 932° 

F give a Prandtl number between 0.006 and 0.007. The 
data points were reevaluated using a Prandtl number of 
0.0065 at 932° F. This increased the Nusselt and Peclet 
numbers of the data by about 16 percent. The reevaluated 
data of Styrikovich and Semenovker are shown in figure 8; 
also shown for comparison is equation (1). { 
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Gilliland, Musser, and Page. — Gilliland, Musser, and Page 
(refs. 2 and 3) measured both heating and cooling coefficients 
for mercury. The heating test section had a 0.319-inch inside 
diameter and a 14-inch length; heat was added by dropwise 
condensation of steam on the outside of the test section. 
The cooling test section had a 0.319-inch inside diameter 
and a 51 -inch length; it was cooled by water flowing on the 
outside in a direction opposite to that of the inside flow. 
Both test sections were made of nickel. The mercury and 
water bulk temperatures entering and leaving the test sec- 
tions and the stream temperature and pressure entering the 
test section were measured. The velocity profile of the mer- 
cury was fully developed at the entrance to both test sections. 
The methods of heating and cooling the mercury were such 
that the heating tests approximated a constant wall tempera- 
ture, while the cooling tests were somewhere between a 
constant wall temperature and a constant heat input. The 
heat-transfer coefficient measured was an over-all average 
coefficient. 

Inasmuch as no wall temperatures were measured, it was 
necessary to separate the mercury heat-transfer coefficients 
from those of the steam and water. This was done by the 
Wilson plot method (see refs. 2 and 3). 

(1) Heating: Tests were run with water in place of mer- 
cury, and the Wilson plot method was used to determine the 
combined resistance of the steam film and the wall. The 
range of water flows covered was sufficiently small and the 
scatter of the points sufficiently great that values of the com- 
bined resistance could be chosen ranging from 40 percent 
greater to 15 percent smaller than the value selected. An 
increase of 40 percent in steam and wall resistance, however, 
would increase the mercury coefficient only about 8 percent. 
An attempt was made to use the results of the mercury runs 
to confirm the steam and wall resistance, but in this case the 
range of mercury flows and data scatter permits selecting a 
value of resistance ranging from 200 percent greater to 50 
percent lower than the value chosen. The slope of the Wilson 
plot for the runs with water can be compared with the slope 
predicted by the standard empirical relation for heat transfer 
to water (ref. 39, p. 168) 


Nu=0.023Re° *Pr aA 


( 8 ) 


The slope predicted by equation (8) turns out to be consider- 
ably higher than the slope best representing the experimental 
data. 

(2) Cooling: At a given mercury flow rate, the water flow 
rate was varied and the combined mercury film and tube wall 
resistance determined by means of a Wilson Plot. The range 
and scatter of the data are such that the resistance of the 
mercury and the wall could be chosen 20 percent lower or 15 
percent higher than the value actually chosen. The corre- 
sponding variation in mercury coefficient would be somewhat 
greater. Alternatively, cooling coefficients for mercury were 
calculated by evaluating the coefficients for water in an 
annulus using the following equation (ref. 39, p. 202): 


StPr 2/3 : 


0 . 020 (DJDi) 0 - 53 
Re °- 2 


( 9 ) 


The resulting mercury coefficients were approximately 40 per- 
cent lower than those derived by the Wilson plot method. 
The physical properties used by Gilliland, Musser, and Page 
are about the same as those of reference 28. In view of the 
possible inaccuracies in the method of evaluating the data, 
the reported results of Gilliland, Musser, and Page may not 
be very accurate. Their data are shown in figure 9 without 
change; shown for comparison are equations (1) and (2). 
The lower values for cooling coefficient may be due to the 
longer length-diameter ratio of the cooling section. 

Elser. — Elser (ref. 4) measured cooling heat-transfer 
coefficients for mercury. Three different test sections were 
used : The test-section inner diameters were 0.317, 0.308, and 
0.260 inch; the 0.317-inch-diameter test section was made of 
mild steel, and the other test sections were made of stainless 
steel. The test sections were all over 38 inches long, but 
measurements were made between two stations 10.2 and 
38.3 inches from the entrance. The mercury was cooled by 
water flowing in a concentric annulus in a direction opposite 
to the flow of mercury. Two thermocouples imbedded in 
the wall measured the wall temperature at the two stations. 
Two other thermocouples immersed in the stream measured 
a temperature close to the fluid centerline temperature. The 
velocity profile of the mercury at the first station was fully 
developed. The cooling-water flow rate was such that a uni- 
form heat input to the wall was approximated. The fully 
developed heat-transfer coefficient was measured. 

The oidy mercury property listed by Elser is Prandtl num- 
ber. These values are in agreement with the values of refer- 
ence 28, and the other mercury properties will be assumed to 
be correct. The basic data are not presented by Elser. He 
presents a plot of Stanton number against Reynolds number 
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Figure 9. — Data of Gilliland, Musser, and Page (refs. 2 and 3) for 
average heat transfer to mercury in round tubes. Length-diameter 
ratio l/D: heating section, 45; cooling section, 160. 


3S5409 — 57 2 




8 


REPORT 1270 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


showing his data points. The values of Stanton number 
have been corrected by Elser to a common Prandtl number 
by approximating the data with a curve of the form 

Nu f —aRe f m Pr f n (10) 

He gives no values of n, so that it is impossible to return to 
the basic points. 

Elser measured mercury flow by measuring the mercury 
pressure drop and assuming the following formula for friction 
factor : 

4/=|p# <80,000 (11) 

4= 2 log (/fe-v'4?)— 0.8 7f<?>80,000 (12) 

v 4 f 

Equation (11) is from Blasius; equation (12) from Karman. 

Elser’s heat-transfer coefficients are based on the difference 
between wall and fluid centerline temperatures. He is not 
certain of the location (depth) of his wall thermocouples and 
states that the difference between a midwall and a wall sur- 
face location results in shifts of heat-transfer coefficients of 
4, 7, and 18 percent, respectively, for the three tubes of 
0.31 7-, 0.308-, and 0.260-inch diameter. In Elser’s data, 
the wall thermocouple is assumed to be at the wall midpoint. 
Martinelli’s predictions for the ratio of the temperature dif- 
ferences (t w —tm)/(t w —tc) (fig. 5(a)) were used to change the 
heat-transfer coefficients of Elser so that they would be 
based on the difference between wall and fluid bulk tempera- 
tures. This increased the Nusselt number about 40 to 60 
percent. In this reevaluation the wall thermocouples were 
assumed to be located at the wall midpoint. If the ther- 
mocouples were assumed at the wall surface, the Nusselt 
numbers would be somewhat increased. Figure 10 shows 



Figure 10. — Reevaluated data of Elser (ref. 4) for fully developed heat 
transfer to mercury in round tubes. 


the reevaluation data of Elser; shown for comparison is 
equation (1). 

Bailey, Cope, and Watson.— Bailey, Cope, and Watson 
(ref. 5) measured cooling coefficients for mercury. The test 
section was a mild-steel tube of 0.437-inch inner diameter. 
The central 18 inches of the tube was surrounded by a water 
jacket, with about 6 inches projecting at each end. These 
ends were enclosed in chambers in such a manner that the 
inlet and outlet mercury passed along the outside of the 
ends before entering and after leaving the test section. Fluid 
temperatures were measured at the inlet and outlet of the 
test section; wall temperatures were measured at four sta- 
tions along the length of the water-jacketed section of the 
tube. 

There is considerable question as to just what temperature 
was measured at the test-section outlet. First, there was 
no provision made for mixing before the exit temperature 
was measured. Second, inasmuch as the mercury was being- 
cooled, the temperature distribution of the mercury was such 
that the temperature near the wall was lower than the bulk 
temperature. The mercury was discharged from the test 
section into a larger chamber, turned 180°, and passed over 
the end of the test section which projected from the water 
jacket. Because of the mixing in the discharge and turning 
processes, the mercury on the outside of the projecting end 
had a nearly flat temperature profile. Hence, the mercury 
on the outside of the projecting end of the test section was 
at about fluid bulk temperature, while the mercury on the 
inside of the projecting end (close to the wall) was at a tem- 
perature lower than fluid bulk temperature. Heat was there- 
fore transferred from the outside to the inside; this tended 
to increase the measured mercury exit temperature and con- 
sequently decrease the observed heat-transfer coefficients. 
The combined effect on heat-transfer coefficient of the heat 
transferred through the projecting end and the lack of mixing- 
before the exit temperature measurement is very difficult to 
estimate. 

The velocity profile at the entrance to the water-jacketed 
section of the test section was close to fully developed. The 
method of cooling was such that uniform wall temperature 
was approximated at the lower mercury Peclet numbers, 
while uniform heat input to the wall was more nearly the 
case at high mercury Peclet numbers. Fully developed heat- 
transfer coefficients were measured. 

The physical properties used by Bailey, Cope, and Wat- 
son were somewhat different from the values of reference 28; 
the Prandtl numbers were about 10 percent high. The data 
of Bailey, Cope, and Watson were therefore reevaluated, 
using the physical properties of reference 28, in two ways: 
First, it was assumed that the measured mercury exit tem- 
perature was equal to the fluid bulk temperature; second, it 
was assumed that the measured mercury exit temperature 
was equal to the fluid centerline temperature, and Martin- 
elli’s predictions of ( t w — (fig. 5) were used to 
calculate the fluid bulk temperature. The results of both 
methods of computation are shown in figure 1 1 ; equations 
(1) and (2) are shown for comparison. Because of the un- 
certainties described in the measurement of mercury exit 
temperature, it. is difficult to say whether either set of data 
in figure 11 is at all correct. 
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Figure 11. — Reevaluated data of Bailey, Cope, and Watson (ref. 5) 
for fully developed heat transfer to mercury in round tubes. 


Lyon. — Lyon (ref. 6) used a tube and concentric annulus 
to measure the combined coefficient resulting from transfer- 
ring heat from a soclium-potassium alloy (52 percent Na, 48 
percent K) flowing in the annulus to the same fluid flowing 
in the tube. The weight flows in the tube and annulus are, 
necessarily, the same. This type of test section is often re- 
ferred to as a “figure eight” and will be so referred to herein. 
Lyon used four different test sections made of nickel and 
having the following dimensions: 


Test section 

A 

B 

c 

D 

Tube inner diameter, D, in 

0.432 

0. 703 

0. 434 

0. 434 

Annulus inner diameter, /?,, in. 

.500 

.757 

.500 

. 500 

Annulus outer diameter, D 0 , in 

.715 

.931 

.684 

.684 

Length, l, in. 

48 

09 

33 

69 


Bulk fluid temperatures were measured at the inlet and out- 
let of the tube and annulus. The velocity profiles of the 
fluid entering the tube and the annulus were approximately 
flat (uniform velocity). The figure-eight test section with 
counterflow gives approximately constant heat input to the 
wall. The heat-transfer coefficients measured were over-all 
average coefficients. 

Lyon used physical properties which were somewhat dif- 
ferent from those of reference 28. The specific heat was 
about 12 percent higher and the thermal conductivity was 
about 6 percent higher. Use of the properties of reference 
28 decreases both the Nusselt and Peclet numbers about 
5 percent. Lyon assumed that the resistances of the walls 
of the four test sections were approximately constant, neg- 
lecting the differences in wall thickness. Lyon did not 
separate the experimental tube and annulus coefficients, but 
rather calculated a combined predicted coefficient using 
equation (1) for the tube and an equation approximating the 
results of Harrison and Menke (ref. 40) 

2Vtt /tM =4.9+0.0175 Pe uan 0.8 (13) 

for the annulus. The Liquid-Metals Handbook (ref. 28) 
mentions equation (13), but prefers equations (4) and (5) 
for heat transfer in an annulus. 

Lyon’s data were reevaluated using the physical properties 
of reference 28 and calculating exactly the resistance of the 


wall. The over-all heat-transfer coefficient was divided into 
a tube coefficient and an annulus coefficient assuming that 
the Nusselt numbers in the tube and annulus are related as 
in equations (6) and (7), which are taken from equations (1), 
(4), and (5). The use of equations (4) and (5) rather 'than 
equation (13) for the annulus results in higher annulus heat- 
transfer coefficients and lower tube heat-transfer coefficients 
for the same over-all heat-transfer coefficient. The re- 
evaluated data of Lyon are shown in figures 12 and 13; 
equations (1), (4), and (5) are shown for comparison. 

TJntermeyer. — The data of Untermeyer were obtained 
from unclassified material in a classified report. Unter- 
meyer measured heating coefficients for a lead-bismuth 
eutectic with and without magnesium addition. The test 
section was a steel tube with a 0.25-inch inner diameter and 
18-inch length. The test section was heated by passing 
electric current directly through it and the fluid it contained. 
Wall temperatures and fluid inlet and outlet temperatures 
were measured. The velocity profile at the test-section 
entrance was closer to flat than to fully developed. The 
method of heating most nearly approximated uniform heat 
input to the wall. Local fully developed coefficients were 
measured. 

The physical properties used by Untermeyer are different 
from those of reference 28. The thermal conductivity used 
by Untermeyer was about 15 percent lower and the vol- 
umetric specific heat was about 8 percent higher. It is 
difficult to determine from the data whether the heat gener- 
ated directly in the fluid has been subtracted from the total 
heat input. It is also difficult to determine whether a mix- 
ing chamber was used in the measurement of the fluid bulk 
temperature leaving the test section. Figure 14 shows the 
data of Untermeyer reevaluated using physical properties 
from reference 28. 



Figure 14. — Reevaluated data of Untermeyer for fully developed heat 
transfer to lead-bismuth eutectic, with and without magnesium 
addition, in round tubes. 
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Figure 12. — Reevaluated data of Lyon (ref. 0) for average heat transfer to sodium-potassium alloy in round tubes. 




Nusselt number, Nu £, Nusselt number, Nu 


REVIEW OF EXPERIMENTAL INVESTIGATIONS OF LIQUID-METAL HEAT TRANSFER 


11 



Pecle! number, Pe 











" 















































































































































































































— 

' 





f ,an~ ^ 

„ 






IF 

(eq. (•< 




-s-" 







Nu 

O TUU£U 

re f.an 

+ )) 




8“ 













— 

— 

— 

— 

— 














o 



Q_ 





<$5 

9 

















1#* 
















c 

O 


y 



















































(b) 





















10 I0 2 I0 3 

Nusselt number, Pe 


Test section A (ratio of outer to inner diameter, 1.43; length- 
diameter ratio l/D, 223). 


(b) Test section B (ratio of outer to inner diameter, 1.23; length- 
diameter ratio l/D, 397). 
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(c) Test section C (ratio of outer to inner diameter, 1.37; length- 
diameter ratio l/D, 179). 


(d) Test section D (ratio of outer to inner diameter, 1.37; length- 
diameter ratio l/D, 375). 


Figure 13.— Reevaluated data of Lyon (ref. 6) for average heat transfer to sodium-potassium alloy in annuli. 
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In round tubes; test section A (length-diameter ratio l/D, 50). 



(b) In annuli; test section A (length-diameter ratio l/D, 54; ratio of 
outer to inner diameter, 1.83). 




(c) In round tubes; test section B (length-diameter ratio l/D, 48). (d) In annuli; test section B (length-diameter ratio l/D, 54; ratio of 

outer to inner diameter, 1.83). 

Figure 15.— Reevaluated data of Werner, King, and Tidball (ref. 7) for average heat transfer to sodium-potassium alloy. 
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Werner, King, and Tidball. — Werner, King, and Tidball 
(ref. 7 and unclassified data from a classified report) used 
a figure-eight test section (tube and concentric annulus 
with same fluid in both) to measure heat-transfer coefficients 
for a sodium-potassium alloy. Cooling coefficients were 
measured in the tube and heating coefficients in the annulus. 
Two test sections having the following characteristics were 
used : 


Test section 

A 

B 


0.68 

0.70 

.75 

1.37 

33.8 

Nickel 




1.37 


33.8 






The tests in test section A were all run with an alloy of 56 
percent sodium and 44 percent potassium. The tests in 
test section B were run with alloys of both 56 percent sodium 
plus 44 percent potassium and 23 percent sodium plus 77 
percent potassium. Fluid temperatures were measured 
at the inlet and outlet of the tube and of the annulus. In 
test section A no provision was made for mixing the fluid 
before measuring the outlet temperatures of the tube or 
the annulus, except that the fluid turned one right-angle 
bend before eacli thermocouple. The outlet temperatures 
measured in test section A were, therefore, somewhere 
between fluid bulk temperature and fluid centerline temper- 
ature, probably closer to fluid centerline temperature. In 
test section B, mixing baffles were used to mix tire fluid 
before measuring outlet temperatures, and the temperatures 
measured were fluid bulk temperatures. The velocity 
profiles of the fluid entering the tube and the annulus were 
essentially flat in test section B, and between flat and fully 
developed in test section A. The figure-eight test section 
with eounterflow gives approximately uniform heat input 
to the wall. The heat-transfer coefficients measured were 
over-all average coefficients. 

Werner, King, and Tidball used physical properties which 
were about the same as those of reference 28. However, 
the relation used to divide the over-all heat-transfer coeffi- 
cient in the test section into separate coefficients for the 
tube and annulus is somewhat different from that recom- 
mended by the Liquid-Metals Handbook (eq. (7)). 

The experimental data of Werner, King, and Tidball were 
reevaluated using equation (7) to separate the over-all 
heat-transfer coefficient into tube and annulus coefficients. 
In addition, the predictions of Martinelli for the ratio of the 
temperature differences {t w —t m ) I (t w —t c ) for the tube and the 
annulus (fig. 5) were used to make allowance for the lack 
of mixing of the fluid before the outlet thermocouples of 
test section A. The reevaluated data of Werner, King, and 
Tidball are shown in figure 15; equations (1) and (5) are 
shown for comparison. 

Sineath.— Sineath (ref. 8) ran heat-transfer tests with 
mercury in rectangular channels. Sineath’s test section 
was of the figure-eight type except that, instead of a tube 
and concentric annulus, he had two rectangular channels 
with one common wall. Heat was added to the mercury in 
one channel and removed from the mercury in the other. 
The common wall of the two channels was 4 inches high 


by y t inch thick and was made of mild steel. The channel 
gap was /( inch and the length, 25 inches. Fluid tempera- 
tures were measured at the inlet and outlet of the two 
channels. No attempt was made to provide any mixing 
of the fluid before the outlet temperatures were measured 
except that the abrupt transition from a 4- by inch 
rectangular channel to the %-inch pipes which carried the 
fluid to and away from the test section probably resulted in 
considerable mixing. The pipe entered the channels at 
right angles to the direction of flow in the channels; there 
was no smooth transition piece between the pipes and the 
channels. The fluid temperatures measured were probably 
close to the bulk temperature. However, the abrupt 
change of section at the entrance to the channels probably 
caused some of the heat-transfer surface to be relatively 
less effective as a result of poor local flow distribution. The 
figure-eight test section with counterflow approximated 
uniform heat input to the wall. The heat-transfer coeffi- 
cients measured were over-all average coefficients. 

Sineath ran four sets of tests. The first three sets were 
inconclusive because of experimental difficulties with air 
entrainment and with deposition of mercurous oxide on 
the wall through which heat was being transferred. These 
problems were partially eliminated in the fourth set of runs. 
There was probably no air entrainment during the fourth 
set of runs; the wall through which heat was being trans- 
ferred was carefully cleaned at the beginning of the runs but 
was covered with a thin layer of scale at the end. 

.Sineath used physical properties similar to those of refer- 
ence 28. The temperatures in the two channels were suffi- 
ciently close that the heat-transfer coefficients in both chan- 
nels could be assumed the same. 

The data of the fourth set of runs of Sineath are shown 
unchanged in figure 16; equation (4) is shown for comparison. 
The data of Sineath are undoubtedly lower than they should 



Peclel number, Pe 


Figure 16. — Data of Sineath (ref. 8) for average heat transfer to mer- 
cury in rectangular ducts. Length-diameter ratio l/D, 50. 
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Figure 17. — Data of English and Barrett (refs. 9 and 10) for full}' 
developed heat transfer to mercury in round tubes. 

be as a result of the deposit of an oxide film on the heat- 
transfer surface and of the abrupt change of cross section 
at the entrance to the channels, which makes a portion of 
the heat-transfer surface ineffective. It is difficult, however, 
to estimate the magnitude of these effects. 

English and Barrett. — English and Barrett (refs. 9 and 10) 
measured heating coefficients for mercury. The test sec- 
tions were of nickel and stainless steel with a 0.051-inch 
inner diameter, a 0.059-inch outer diameter, and a 1.9-inch 
length. A copper coating was bonded to the outside of the 




Figure 18. — Data of English and Barrett (ref. 9) for entrance-region 
heat transfer to mercury in round tubes. Peclet number, 800 to 

900. 


test section; the outer diameter of the copper was 0.0825 
inch. The test section was heated by passing electricity 
directly through it. The inlet and outlet mercury bulk 
temperatures were measured, as was the outside-wall temper- 
ature along the test section; the voltage distribution along 
the test section was also measured. The velocity profile 
at the test-section entrance was fully developed. The 
method of heating most nearly approximated uniform heat 
input to the wall. English and Barrett measured local heat- 
transfer coefficients along the test section and present the 
local fully developed coefficients for all runs. For one run, 
the local coefficient along most of the tube is presented. 

The physical properties used by English and Barrett 
are the same as those of reference 28 except that the viscosity 
is slightly high at low temperatures. This will probably not 
affect the fully’- developed heat-transfer coefficients, but the 
entrance-region Reynolds numbers should be increased 
3 to 4 percent. 

The fully developed heat-transfer coefficients of English 
and Barrett are shown unchanged in figure 17; equation (1) 
is shown for comparison. The entrance heat-transfer 



Figure 19. — Data of Seban (ref. 11) for fully developed heat transfer 
to lead-bismuth eutectic in round tubes. 


Figure 20. — Data of Seban (ref. 11) for entrance-region heat transfer 
to lead-bismuth eutectic in round tubes. 
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coefficients for the one run presented are shown in figure 
18; Deissler’s predicted curves for a Prandtl number of 0.01 
and the same Peclet number range are shown for comparison. 

Seban. — Seban (ref. 11) measured heat-transfer coefficients 
with lead-bismuth eutectic in two different types of test 
section. One was the figure-eight type with a tube and 
concentric annulus; the other was a copper-coated tube 
heated by external electric heaters. Only combined heat- 
transfer coefficients for the tube and annulus of the figure- 
eight test section are presented in reference 11. Not enough 
basic data (specifically, fluid temperatures) are presented to 
separate the tube and annulus coefficients. Accordingly, 
only the electrically heated test section will be discussed. 
The test section had a 0.652-inch inner diameter and a 48- 
inch length. The copper coating was for the purpose, of 
containing the wall thermocouples in a region of relatively 
low temperature gradient and of smoothing out the non- 
uniformities of heat input of the external electric heaters. 
The fluid bulk temperatures were measured at the inlet 
and outlet of the test section, and the wall temperatures 
were measured at eight stations along the tube. The velocity 
profile was close to fully developed at the entrance to the 
test section. The method of heating approximates very 
closely uniform heat input to the wall. Local heat-transfer 
coefficients were measured. The local fully developed 
coefficients are presented for all the runs; entrance coeffi- 
cients are presented for a few of the runs. 

The physical properties of Seban are the same as those 
given in reference 28. Seban had some trouble with fouling, 
which caused the heat-transfer coefficients to decrease with 
time. Figure 19 shows unchanged the fully developed heat- 
transfer coefficients of Seban; equation (1) is shown for 
comparison. Those points taken immediately after cleaning 
have higher heat-transfer coefficients than the others. 
The entrance heat-transfer coefficients presented by Seban 
are shown in figure 20; predicted curves of Deissler (see 
fig. 8) for the same range of Peclet number and for a Prandtl 
number of 0.01 are shown for comparison. 

Trefethen. — Trefethen (refs. 12 and 13) used a figure-eight 
type of test section to measure heat-transfer coefficients with 
mercury. Six different tubes, described in the following 
table, were used in the tube and concentric annulus test 
section: 


Test section 

A 

B 

C 

D 

E 

F 

Tube inner diam- 

0.711 

0.737 

0.585. . . 

0.523. .. . 

0.308 . . 

0.429. 

eter, D, in. 







Annulus inner di- 

.749. . . 

.748-..- 

.629. . 

.627-. . 

.378 .. 

.500. 

ameter, Di, in. 







Annulus outer di- 

.874 

.874 

.874- - . 

.874 

.874 -- 

.874. 

ameter, Z> 0 , in: 







Length, l, in. ... 
Tube material-.. 

39.3 

39.3 

39.3 

39.3 

39.3. .. 

39.3. 

Stainless 

steel. 

Stainless 

steel. 

Stainless 

steel. 

Copper. 

Copper.. 

Copper. 


Heating and cooling tests were run in both the tube and 
the annulus. Trefethen measured the fluid bulk tempera- 
ture at the inlet and outlet of the tube and annulus. He 
also measured the wall temperature of the outside of the 
annulus. The velocity profiles were between flat and fully 
developed at the entrance to the test section, probably a 
little closer to flat. The counterflow figure-eight test 
section approximated a uniform heat input to the wall. 

385409. — 56- 3 

X 


Trefethen presents a fully developed heat-transfer coef- 
ficient for the central section of his tube (from 10.2 to 29.4 
in. from the entrance.) In his calculations he assumed that, 
the temperature difference between the fluid in the tube 
and in the annulus remains the same as the temperature 
difference at the entrance to the tube and annulus, and that 
the fluid bulk temperature gradient along the length of the 
tube center section is the same as the temperature gradient 
along the annulus outer-wall center section. Trefethen 
separates the tube and annulus coefficients in a manner 
different from that resulting from the use of equations (6) 
and (7). 

The physical properties used by Trefethen are about the 
same as those of reference 28 for the range of temperature 
covered by his experiments. (The values of thermal con- 
ductivity at high temperature (extrapolated by Trefethen 
to correct the data of Styrikovich and Semenovker) are 
lower than those of the Liquid-Metals Handbook by about 
6 percent at 212°, 14 percent at 392°, and 20 percent at 
662° F.) 

The data of Trefethen for fully developed heat-transfer 
coefficients were recalculated using equations (6) and (7) 
to separate the coefficients of the tube and annulus. The 
reevaluated data for the tube are shown in figure 21 ; equation 
(1) is shown for comparison. Trefethen gives enough data 
to permit the calculation of over-all average heat-transfer 
coefficients for the tube and annulus. These coefficients 
also were calculated using equations (6) and (7) to separate 
the individual coefficients in the tube and annulus. These 
data are shown in figure 22; equations (1), (4), and (5) 
are shown for comparison. 

Doody and Younger. — Doodv and Younger (refs. 14 and 
15) measured both heating and cooling coefficients for 
mercury with and without sodium additions. The test 
section was a steel tube 0.493 inch in inner diameter and 61 






























































Test inner diameter, 

section D 0 /Di 

A, B, C, D <1.4 
























o 












O 

F 


1 

7 

5 



























7 







Nu f 

= 7.0 + O.C 

)25 

Pi 

V 

xi 

Hi 


.(!» _ 
















- 

- 

- 

-- 

i 


8® 

0 












— ■ 




□ 

s 











- 

■ 





£ 

















oM 


r 

















\ ■ ■ 














“1 






ro 

















































































































J 






i i 1 i i i i i i i i i i i i i i 

I0 2 I0 3 I0 4 


Peclet number, Pe 


Figure 21. — Reevaluated data of Trefethen (refs. 12 and 13) for fully 
developed heat transfer to mercury in round tubes. 
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(a) In round tubes; length-diameter ratio l/D, 53 to 128. 
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(b) In annuli; length-diameter ratio l/D, 79 to 334. 

Figure 22. — Reevaluated data of Trefethen (refs. 12 
and 13) for average heat transfer to mercury. 


inches long. The test section was heated or cooled by water 
flowing in a concentric annulus. Both parallel and counter- 
flow runs were made. Theannulus was 61 inches long, but 
the annulus entrance and exit were each 6 inches from the 
ends of the test section; therefore, the length of the test 
section between the annulus entrance and exit was 49 
inches. Tube wall temperatures were measured at five 
stations starting % inch downstream of the annulus inlet 


and ending 3 inches upstream of the annulus outlet (inlet 
and outlet refer to the parallel-flow case) ; the measurements 
covered 45% inches of the test seection. The mercury 
temperature was measured at the inlet and outlet of the 
test section. The measurement of test-section exit tem- 
perature was made without any preliminary mixing of the 
fluid, and the exit temperature measured is closer to the 
fluid centerline temperature than to the fluid bulk tempera- 
ture. The mercury velocity profile at the test-section 
entrance was probably closer to flat than it was to fully 
developed. The method of heating resulted in a wall 
condition somewhere between uniform heat input and uni- 
from wall temperature for the counterflow runs. The 
parallel-flow runs resulted in a wall condition where the 
rate of heat input varied even more rapidly than for the 
condition of uniform wall temperature. 

The physical properties used by Doody and Younger are 
the same as those of reference 28 except for the thermal 
conductivity. The values of thermal conductivity used by 
Doody and Younger are low by about 1 to 14 percent in the 
temperature range of the investigation. Because of the 
location of the annulus entrance and exit as described, some 
effective length of test section between 49 and 61 inches must 
be selected. Doody and Younger used a method due to 
Sherwood and Petrie (ref. 41) and arrived at an effective 
length of 56 inches. Since the wall temperatures at the ends 
of the test section were not measured, Doody and Younger 
extrapolated the wall temperature measurements to cover a 
length of 56 inches, the “effective” length of their test sec- 
tion. Inasmuch as this extrapolation necessarily neglects 
end effects, the data of Doody and Younger represent some- 
thing between over-all average and fully developed heat- 
transfer coefficients. 

The heat balances of Doody and Younger show deviations 
as great as 140 percent, with deviations between 40 and 100 
percent being quite common. Flow was measured by an 
orifice, and the orifice calibration showed variations as great 
as 50 percent. The end temperature differences between 
the wall and the fluid found by the previously mentioned 
extrapolation were very small, varying from about 0.3° to 
8° F, with values of 2° F or less being extremely common. 
Small errors in temperature measurement can therefore 
result in large errors in log mean temperature difference. 

Doody and Younger attempted to check their experimental 
apparatus by running heat-transfer experiments with butanol. 
Unfortunately, most of these data were in the transition 
region. Some of the data were in the laminar-flow region, 
and these data were 25 to 75 percent higher than the predic- 
tions of the Colburn equation for laminar flow (ref. 39, p. 191). 

In view of the difficulties mentioned, the data of Doody 
and Younger may not be very accurate. The data of 
Doody and Younger were reevaluated using the physical 
properties of reference 28 and the predictions of Martinelli 
for the ratio (#«, — t m )/(t w — t c ) to determine the value of the 
temperature differences between the wall temperature and 
the fluid bulk temperature at the test-section exit. The 
original data on wall temperature are not presented in either 
reference 14 or 15; therefore the extrapolated end tern per a- 
i ture will be used for the wall temperature at the ten-section 



REVIEW OF EXPERIMENTAL INVESTIGATIONS OF LIQUID-METAL HEAT TRANSFER 


17 



Figure 23. — Reevaluated data of Doody and Younger (refs. 

inlet, and outlet. The reevaluated data of Doody and 
Younger are shown in figure 23; equations (1) and (2) are 
shown for comparison. 

Lubarsky. — Lubarsky (ref. 16) used a figure-eight type of 
test section to measure heating coefficients in a tube and 
cooling coefficients in a concentric annulus for lcacl-bismuth 
eutectic with and without magnesium additions. The test 
section was 40.2 inches long, with a 0.402-inch tube inner 
diameter, a 0.50-inch annulus diameter, and a 0. 625-inch 
annulus outer diameter. The fluid bulk temperatures at 



Peclet number, Pe 

(b) With small sodium additions. 

14 and 15) for heat transfer to mercury in round tubes. 


the inlet and outlet of the tube and annulus were measured. 
The entering velocity profile was approximately flat. The 
figure-eight-type heat exchanger with counterflow approxi- 
mated uniform heat input to the wall. Over-all average 
heat-transfer coefficients were measured. 

Lubarsky used physical properties which were the same 
as those of reference 28. He used equation (6) to separate 
the heat-transfer coefficients in the tube and the annulus. 
Lubarskv’s data are shown unchanged in figure 24; equations 
(1) and (4) are shown for comparison. 
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(a) In round tubes; length-diameter ratio IjD, 100. 



Figure 24. — Data of Lubarsky (ref. 16) for average heat transfer to lead-bismuth eutectic. 
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Johnson, Hartnett, and Clabaugh (lead-bismuth tests). — 
Johnson, Hartnett, and Clabaugh (refs. 17 and 18) measured 
heating coefficients for lead-bismuth eutectic in an aluminum- 
coated tube heated externally by electric heaters. The test 
section was very similar to that of Seban (ref. 11) described 
previously, except that an aluminum coating was used instead 
of a copper one. The test-section inner diameter was 0.652 
inch and its length 48 inches. The fluid bulk temperature 
was measured at the inlet and outlet of the test section, and 
the wall temperature was measured at eight stations along 



Figure 25. — Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 
18) for fully developed heat transfer to lead-bismuth eutectic in 


round tubes. 



(a) Length-diameter ratio x/D, 4.6. 


the wall. The inlet velocity profile was very close to fully 
developed. The method of heating very closely approxi- 
mated uniform heat input to the wall. Local heat-transfer 
coefficients are presented for both the fully developed region 
and the entrance region. 

Physical properties the same as those of reference 28 were 
used. The data on fully developed heat-transfer coefficients 
are shown unchanged in figure 25; equation (1) is shown for 
comparison. The data on entrance heat-transfer coefficients 
are shown unchanged in figure 26 ; predicted curves of Deissler 
(fig. 3(b)) for the same length-diameter ratio x/D and a 



Peclet number, Pe 


(b) Length-diameter ratio x/D, 13.8. 



(c) Length-diameter ratio x/D, 23. 


Figure 26.— -Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 18) for entrance-region heat transfer to lead-bismuth eutectic in round tubes. 
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Peclet number, Pe 


Figure 27.-- Data of Johnson, Hartnett, and Clabaugh (refs. 17 and 

18) for average heat transfer to lead-bismuth eutectic in round 

tubes. Length-diameter ratio l/D, 74. 

Prandtl number of 0.01 are shown for comparison. The pre- 
dictions of Poppendiek and Palmer (fig. 3(b)) for low Reyn- 
olds numbers are also shown in figure 26. From the local 
data of Johnson, Hartnett, and Clabaugh, it is possible to 
determine an approximate over-all average coefficient by 
plotting and integrating the local coefficients. The results 
of this procedure are shown in figure 27; equation (1) is 
shown for comparison. 

Isakoff and Drew. — Isakoff and Drew (refs. 19 and 20) 
measured heating coefficients with mercury. The test sec- 
tion was a stainless-steel tube with a 0.127-inch wall thickness, 
1 .5-inch inner diameter, and about 223-inch length, heated 
externally by electric heaters. The fluid bulk temperature 
was measured at the test-section inlet and outlet, and the 



outside wall temperature was measured at seven stations 
along the tube. Velocity and temperature profiles in the 
fluid were measured at three stations along the tube (x/JJ= 
58, 98, and 138). The entrance velocity profile was very 
close to flat. The method of heating approximated very 
closely uniform heat input to the wall. Local fully de- 
veloped heat-transfer coefficients were measured at the 
stations of x/D=98 and x/D= 138. The heat-transfer coeffi- 
cients measured at the x/D=5S are still in the entrance 
region. 

The physical properties used by Isakoff and Drew are the 
same as those of reference 28. The inside wall temperature 
was calculated in two ways: one was to extrapolate the 
temperature profile in the fluid to the wall ; the other was to 
use the measured outside wall temperature to calculate the 
temperature drop through the wall. When this calcula- 
tion was made, the two inside wall temperatures were found 
to coincide for only three of the total of 12 experimental runs. 
For the other nine runs, the inside wall temperature cal- 
culated from the outside wall temperatures was higher than 
the inside wall temperature as extrapolated from the fluid 
temperature profile. Fluid and wall temperatures for two 
typical runs are shown in figure 28. 

This discrepancy between the two methods of determining 
inside wall temperature may be due to inaccuracies in the 
measurement of outside wall temperature. The outside wall 
temperature was measured by eight thermocouples at each 
station. The temperature readings of these thermocouples 
varied as much as 25 percent from each other in the high-flux 
region. The deviation may have been due to the proximity 
of the thermocouples to the electric heaters. At any rate, 
the order of magnitude of the variation of the thermocouple 
readings on the outside wall is as great as the magnitude of the 
differences in temperature resulting from the two methods of 
calculating inside wall temperature. It is interesting to note, 
however, that the inside wall temperature as calculated 
from the outside wall temperature is greater (for nine cases 
out of twelve) than that extrapolated from the fluid tempera- 
ture profile. This effect is that which would be noted if 



0 2 4 6 8 10 12 14 16 18 

Distance along test section, ' x , ft 

(b) Run 12. 


Figure 28.— Test-section temperatures from two typical runs of Isakoff and Drew (refs. 19 and 20). 
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results; Deissler’s predicted curves (fig. 3) are shown for 
comparison. The outside wall temperature was measured 
at a sufficient number of stations along the tube (see fig. 28) 
to permit the estimate of an over-all average heat-transfer 
coefficient. The average over-all coefficient from x/D=Q.3 
to x/D= 138, based on inside wall temperatures calculated 
from the outside temperature, is shown in figure 31. These 
average coefficients are actually lower than the fully devel- 
oped coefficients because the average outside wall tempera- 
ture at the second station from the test-section entrance is 
higher than might be expected from the other measured 
temperatures. Whether this might be a result of the local 
temperature gradients caused by the external electric heaters 
cannot be determined. 

The temperature profiles in the fluid at x/D=138 are 
shown in figure 32 and compared with the predictions of 
Martinelli (ref. 27). 

Stromquist. — Stromquist (ref. 21) measured heating co- 
efficients for mercury with and without sodium additions. 
The test section was a steel tube heated by passing electricity 
directly through the tube. The following four different test 
sections were used: 


Figure 29. — Data of Isakoff and Drew (refs. 19 and. 20) for fully 
developed heat transfer to mercury in round tubes. Length- 
diameter ratio a ',/D, 138. 

there were some form of interfacial resistance between the 
fluid and the tube. However, because of the circumferen- 
tial variation of outside wall temperature and because three 
of the runs showed no difference in the insidewall temperature 
calculated by the two methods, no conclusions can be reached. 

The data of Isakoff and Drew are shown in figure 29 for 
the fully developed heat-transfer coefficient (x/D= 138); the 
coefficient is shown for both methods of calculating inside 
wall temperature; equation (1) is shown for comparison. 
The entrance data at x/D = 58 are shown in figure 30; in- 
asmuch as the ordinate in this figure is a ratio, both methods 
of calculating the inside wall temperature give the same 



Figure 30. — Data of Isakoff and Drew (refs. 19 and 20) for entrance- 
region heat transfer to mercury in round tubes (inside wall temper- 
atures extrapolated from fluid temperature profile). Length- 
diameter ratio xjD , 58. 


Test section 

A 

B 

C 

D 

Inner diameter, D { , in. 

0. 380 

0. 488 

0. 753 

0. 787 

Outer diameter, D 0 , in. 

. 753 

1.002 

1.500 

1. 501 

Length, l, in 

47. 25 

48.25 

50. 25 

50. 25 


Fluid bulk temperatures were measured at the test-section 
inlet and outlet, and the outside wall temperatures were 
measured at 12 stations along the length of the test section. 



Figure 31. — Data of Isakoff and Drew (refs. 19 and 20) for average 
heat transfer to mercury in round tubes (inside wall temperatures 
calculated from outside wall temperatures). 
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Figure 32. — Data of Isakoff and Drew (refs. 19 and 20) for fully 
developed temperature distribution for heat transfer to mercury 
in round tubes. Length-diameter ratio x/D, 138. 

The entrance velocity profile to the test section was some- 
where between flat and fully developed. The method of 
heating more nearly approximated uniform heat input than 
it did uniform wall temperature. Local heat-transfer co- 
efficients, both entrance and fully developed, were measured. 

Stromquist used physical properties which were the same 
as those of reference 28. Figure 33 shows the fully developed 
heat-transfer coefficients of Stromquist unchanged; equation 
(1) is included for comparison. The unchanged entrance 
heat-transfer coefficient data of Stromquist are shown in 
figure 34; Deissler’s predicted curves (fig. 3(a)) are shown 
for comparison. The predictions of Poppendiek and Palmer 
(fig. 3(a)) for low Reynolds numbers are shown in figure 
34(a). 

MacDonald and Quittenton.— MacDonald and Quittenton 
(refs. 22 and 23) measured heating coefficients with sodium. 
The test section consisted of a monel tube with a copper 
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Figure 33. — Data of Stromquist (ref. 21) for fully developed heat 
transfer to mercury, with and without sodium additions, in round 
tubes. 

jacket bonded to the outside. The test section of 0.625- 
inch inner diameter and 60.05-inch length was heated ex- 
ternally by electric heaters. The purpose of the copper 
jacket was the same as the purpose of the copper coating- 
used by Seban (ref. 11) and described in a previous section 
entitled “Seban.” The fluid bulk temperature was measured 
at the inlet and outlet of the test section, and the wall tem- 
peratures were measured at 1 1 stations along the test sec- 
tion. The entrance velocity profile was close to fully de- 
veloped if the length of piping in the diagram shown in 
figure 1 of reference 23 is drawn to scale. The method of 
heating very closely approximated uniform heat input to 
the walls. Local heat-transfer coefficients were measured; 
the authors present fully developed heat-transfer coefficients 
for a length of the tube from 47.3 to 54.8 inches down- 
stream of the tube entrance. 

MacDonald and Quittenton used the same physical proper- 
ties as those in reference 28. The data for fully developed 
heat-transfer coefficients are shown unchanged in figure 
35; equation (1) is shown for comparison. The data show 
a great amount of scatter. Consecutive runs at identical 
Peclet numbers and similar temperature levels vary as much 
as 120 percent in Nusselt number, with variations of 30 to 
60 percent in consecutive runs being common. In view of 
this scatter, the entrance coefficients and over-all average 
coefficients have not been calculated, although the data were 
sufficient to make these calculations possible. 

Johnson, Clabaugh, and Hartnett (mercury tests).— 
Johnson, Clabaugh, and Hartnett (ref. 24) measured heat- 
ing coefficients for mercury. The test section was almost 
identical to the test section described previously under the 
section “Johnson, Hartnett, and Clabaugh (lead-bismuth 
tests).” The test section was an aluminum-coated tube 
with a 0.652-inch inner diameter and a 48-inch length. 
The fluid bulk temperatures were measured at the inlet and 
outlet of the test section; the wall temperatures were meas- 
ured at eight stations along the test section. The entrance 
velocity profile was close to fully developed. The method of 
heating approximated very closely uniform heat input to 
the wall. Local heat-transfer coefficients, both fully devel- 
oped and entrance, are presented. 
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The physical properties used by Johnson, Clabaugh, and 
Hartnett are the same as those of reference 28. The fully 
developed heat-transfer coefficients of Johnson, Clabaugh, 
and Hartnett are shown unchanged in figure 36; equation 
(1) is shown for comparison. The entrance heat-transfer 
coefficients are shown in figure 37. Deissler’s predicted 
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Figure 35. — Data of MacDonald and Quittenton (refs. 22 and 23) 
for fully developed heat transfer to sodium in round tubes. 




(e) 5000<Peclet number<20,000. 


Figure 34.— Data of Stromquist (ref. 21) for entrance-region heat transfer to mercury, with and without sodium additions, in round tubes. 
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Peclet number, Pe 

Figure 36. — Data of Johnson, Clabaugh, and Hartnett (ref. 24) for 
fully developed heat transfer to mercury in round tubes. 

curves (fig. 3(b)) arc shown for comparison; the predictions 
of Poppendiek and Palmer (fig. 3(b)) for low Reynolds 
number are also shown in figure 37. From the experimen- 
tally determined entrance and fully developed heat-transfer 
coefficients, it is possible to determine by integration the 
over-all average coefficient. The resulting over-all average 
heat-transfer coefficients are shown in figure 38; equation 
(1) is shown for comparison. 



(a) Length-diameter ratio x/D, 4.6. 


Johnson, Hartnett, and Clabaugh (laminar and transition 

flow). — Johnson, Hartnett, and Clabaugh (ref. 25) have 
measured heating coefficients for lead-bismuth eutectic and 
mercury in the laminar and transition flow regions. The 
test section used was identical to the test sections used in 
the investigations of lead-bismuth eutectic and mercury in 
the turbulent-flow region by Johnson, Hartnett, and Cla- 
baugh (see the preceding sections entitled “Johnson, Hart- 
nett, and Clabaugh (lead-bismuth tests)” and “Johnson, 
Clabaugh, and Hartnett (mercury tests)”). The test sec- 



(b) Length-diameter ratio x/D, 13.8. 



Figure 37. — Data of Johnson, Clabaugh, and Hartnett (ref. 24) for entrance-region heat transfer to mercury in round tubes. 
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Figure 38. — Data of Johnson, Clabaugh, and Hartnett (ref. 24) for 
average heat transfer to mercury in round tubes. Length-diameter 
ratio l/D, 74. 


tion was an aluminum-coatecl tube of 0.652-inch inner 
diameter and a 48-inch length. The fluid bulk tempera- 
tures were measured at the inlet and outlet of the test sec- 
tion, and the wall temperatures were measured at eight sta- 
tions along the test section. The entrance velocity profile 
was in doubt, since the flow was mostly in the transition 
region. The method of heating approximated very closely 
uniform heat input to the wall. Local heat-transfer co- 



Figure 40. — Data of Poppendiek and Harrison (ref. 26) for average 
heat transfer to mercury in very short round tubes. 


efficients, both fully developed and entrance, were measured. 

The physical properties used by Johnson, Hartnett, and 
Clabaugh are the same as those of reference 28. The fully 
developed heat-transfer coefficients are shown in figure 
39(a) ; equation (1) is shown for comparison. The entrance- 
region heat-transfer coefficients are shown in figure 39(b). 
There has been no theoretical work on entrance-region heat 
transfer in the transition flow region; therefore no curves 
can be shown for purposes of comparison. 
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(a) Fully developed heat transfer. (b) Entrance-region heat transfer. 

Figure 39. — Data of Johnson, Hartnett, and Clabaugh (ref. 25) for heat transfer to lead-bismuth eutectic and mercury in round tubes in 

laminar and transition flow regions. 
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Poppendiek and Harrison. — Poppendiek and Harrison 
(ref. 26) have measured average heating coefficients with 
mercury in very short test sections. The test section was a 
small hole along the axis of a copper disk of 3-inch outer 
diameter heated on the outside with water. An unheated 
starting length was used so that the entrance velocity profile 
was very close to fully developed. Three different test sec- 
tions were used: 


Test section 

A 

B 

C 

Inner diameter, Di, in --- 

h 

Mo 


Length, l, in 

w 

Me 

Mo 



(a) 600<Peclet number<1000. 



(c) 2000< Peelet number<3000. 


The fluid bulk temperatures of the mercury were measured 
at the inlet and outlet of the test section. Wall tempera- 
tures were measured at several radial stations in the test 
section. The method of heating approximated constant wall 
temperature. 

The properties used by Poppendiek and Harrison are the 
same as those of reference 28. The over-all average heat- 
transfer coefficients are shown in figure 40. Also shown is 
a predicted curve of Poppendiek and Harrison (ref. 26) for 
average coefficients. They obtained the curve by integrat- 
ing the local coefficients predicted by Poppendiek and 
Palmer for low Peelet numbers (fig. 3). 



(b) 1000< Peelet number<2000. 



Figure 41. — Data of Poppendiek and Harrison (ref. 26) for average heat transfer to mercury in very short round tubes compared with 

predictions of Deissler (ref. 37). 




Nusselt number, Nu 
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Isokoff and Drew (mercury; inside wall tempera- 
tures calculated from fluid temperature 
profiles), refs 19 and 20 
Isakoff and Drew (mercury; inside wall tempera- 
calculated from outside wall temperatures), 
refs. 19 and 20 

Johnson, Hartnett, and Clabaugh (lead-bismuth 
eutectic), refs 17 and 18 
Sty rikovich and Semenovker (mercury), ref. I 
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(a) Faired curves of experiments. 


So that the data could be compared with the predictions 
of Deissler (fig. 3), the local heat-transfer predictions of 
Deissler were integrated to obtain predictions for average 
Nusselt number for short length-diameter ratios. In figure 
41 the data of Poppendielc and Harrison for average Nusselt 
number are compared with the predictions of Deissler for 
the average Nusselt number for small length-diameter ratios. 

INTERCOMPARISON OF EXPERIMENTAL RESULTS AND 
COMPARISON WITH THEORY 

The experimental results of the various investigators will 
be compared with each other and with theoretical predic- 
tions. The arrangement of the subjects to be considered 
will be the same as in the section Theoretical Investigations 
of Liquid-Metal Heat Transfer. 

Fully developed heat-transfer coefficients. — Fullv devel- 
oped heat-transfer coefficients for the case of uniform heat 
input to the wall were measured by the following investi- 
gators from the group of 20 investigations reviewed; Sty- 
rikovich and Semenovker; Elser; Untermeyer; English and 
Barrett; Seban; Trefethen; Johnson, Hartnett, and Clabaugh 
(lead-bismuth tests); Isakoff and Drew; Stromquist; Mac- 
Donald and Quittenton; Johnson, Clabaugh, and Hartnett 
(mercury) ; and Johnson, Hartnett, and Clabaugh (laminar 
and transition flow). Curves representing mean lines 
through the data of these various investigators are 
shown in figure 42(a); when the amount of scatter 
of a set of data is so great that no mean line 
can be drawn through it, a cross-hatched area is used 
to represent the data; equation (1) is shown in figure 42(a) 
for purposes of comparison. The spread of all the data in 
figure 42(a) is extremely great. However, the following 
data are not considered; 

(1) Data below a Peclet number of 200: These data, being 


(b) 


- Lyon (theo.) , ref. 6 (eg. (I) ) 

Trefethen (mercury), refs. 12 and 13 
Johnson, Hartnett, and Clabaugh (mercury and lead-bismuth eutectic; 

laminar and transition flow), ref. 25 
Johnson, Clabaugh, and Hartnett (mercury), ref. 24 
Stromguist (mercury), ref. 21 
English and Barrett (mercury) refs. 9 and 10 
Untermeyer (lead-bismuth eutectic) 

Untermeyer ( lead-bismuth eutectic plus magnesium) 

Seban (lead-bismuth eutectic), ref. II 

Isakoff and Drew (mercury, inside wall temperatures calculated 
from fluid temperature profiles), refs. 19 and 20 
Isakoff and Drew (mercury, inside wall temperatures calculated 
from outside wall temperatures), refs. 19 and 20 
Johnson, Hartnett, and Clabaugh (lead-bismuth eutectic), refs. I7and 18 
Styrikovich and Semenovker (mercury) ref. I 
MacDonald and Quittenton (sodium), refs. 22 and 23 
Elser (mercury), ref. 4 




-6— O- 


kl 






IO 3 

Peclet number Pe 


Nu- 0.625 Pe 0A 


(b) Reevaluated data points. 

Figure 42. — Comparison of measured and predicted fully developed Nusselt numbers in round tubes with constant heat input to wall. 
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below a Reynolds number of 10,000 and therefore in the 
transition flow region, are not intended to be represented by 
equation (1). 

(2) Data of Elser, and MacDonald and Quittenton: These 
are less reliable because of the very large scatter of the 
data. 

(3) Data of Untermeyer: Severe corrosion throughout 
the duration of the tests caused large changes in the physical 
dimensions of the test section, as well as possibly contamin- 
ating the fluid and the heat-transfer surface. 

The Nusselt number of the remaining data can be compared 
with the predicted values of equation (1) as follows: 


Peclet number, Pe 

200 

500 

1000 

2000 

5000 

9000 

Range of ratio of measured values . 

0.54 j 

0. 55 

0.57 : 

0. 57 

0. 55 

0. 50 

to predicted values. 

to 

to 

to 

to 

to 

to 


0.69 

0. 75 

0. 96 

1.04 

1. 13 

1. 26 


Another method of comparing the data for the fully 
developed heat transfer is to show on a single plot the actual 
corrected data of all the investigators (fig. 42(b)). If the 
same data are considered valid in this figure as in figure 
42(a), a line given by the following equation would best 
represent most of the data: 

Nu= 0.625 Pe 0i (14) 

This equation is purely empirical and does not in any way 
suggest that the theoretical predictions are faulty. How- 
ever, inasmuch as there is a considerable amount of scatter 
and since most of the data agree fairly well with this line, 
it would seem preferable for the designer to use equation (14) 
until further experiment reduces the uncertainty as to the 
precise values of liquid-metal heat-transfer coefficients. 

Local heat-transfer coefficients in entrance region. — 
Entrance-region heat-transfer coefficients have been meas- 
ured by the following investigators from the group of 20 
investigations reviewed: English and Barrett; Seban; John- 
son, Hartnett, and Clabaugli (lead-bismuth eutectic); 
Isakoff and Drew; Stromquist; Johnson, Clabaugli, and 
Hartnett (mercury); Johnson, Hartnett, and Clabaugli 
(laminar and transition flow); and Poppendiek and Harrison. 
The bulk of the data on heat transfer in the entrance region 
is in the reports of Johnson, Hartnett, and Clabaugli (lead- 
bismuth eutectic); Stromquist; Johnson, Clabaugli, and 
Hartnett (mercury); and Poppendiek and Harrison; and is 
presented in figures 26, 34, 37, and 41. 

There is considerable scatter in most of the entrance 
heat-transfer data presented. The predictions of Deissler 
(fig. 3) agree well with the data of Stromquist (fig. 34), but 
fall slightly low when compared with the remaining data 
(figs. 26, 37, and 41). 

English and Barrett, Seban, and Isakoff and Drew present 
a small amount of entrance-region heat-transfer data. As 
may be seen from figure 18, the data of English and Barrett 
are considerably lower than the predictions of Diessler. 
The data of Seban (fig. 20) agree reasonably well with the 
predictions of Deissler. The data of Isakoff and Drew 
(fig. 30) are considerably higher than the predictions of 
Deissler, which may be, in part, due to the fact that the 
entrance velocity profile of Isakoff and Drew was very 


nearly flat, while in the analyses of Deissler a fully developed 
velocity profile, was assumed at the entrance. 

Average heat-transfer coefficients. — Theoretical predic- 
tions of the over-all average heat-transfer coefficient can be 
made from the information on local heat-transfer coef- 
ficients. The predictions of Deissler for the ratio Nu r /Nu x 
were integrated mechanically, and the values of the ratio of 
average Nusselt number to fully developed Nusselt number 
NuJNu r are shown plotted against length-diameter ratio 
x/D for various Peclet numbers in figure 43. 

Values of average Nusselt number were determined from 
the values of the ratio Nu av /Nu f in figure 43 and the values 
of fully developed Nusselt number Nu r of equation (1). 
The results are shown in figure 44, which gives the variation 
of average Nusselt number with Peclet number for several 
length-diameter ratios. 



Length- diameter ratio, x/D 


Figure 43. — Predictions of Deissler (ref. 37) for variation of ratio of 
average Nusselt number to fully developed Nusselt number with 
length-diameter ratio for various Peclet numbers (Prandtl number, 
0 . 01 ). 



Peclet number, Pe 


Figure 44. — Variation of average Nusselt number with Peclet number 
for several length-diameter ratios as determined from figure 43 and 
equation (1). 
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Figure 45. — Comparison of measured and predicted average Nusselt 
numbers in round tubes with uniform heat input to wall. 


The measured average heat-transfer coefficients are 
described as follows: 

(1) Uniform heat input to the wall; round tubes: Average 
heat-transfer coefficients in round tubes with constant heat 
input to the walls were measured by Lyon; Werner, King, 
and Tidball; Trefethen; Lubarsky; Johnson, Hartnett, and 
Clabaugh (lead-bismuth); Isakoff and Drew; and Johnson, 
Clabaugh, and Hartnett (mercury). Curves representing 
mean lines through the data of these various investigators 
are shown in figure 45. Also shown is the relation for average 
Nusselt number (l/D= 100) from figure 44. The Nusselt 
numbers of the data compare with the predicted values for 
a length-diameter ratio of 100 as follows (values below a 
Peclet number of 200 are not considered because they 
fall in the transition flow region): 


Peclet number, Pe 

200 

500 

1000 

2000 

5000 

9000 

Range of ratio of measured values 

0.64 

0. 75 

0.54 

0.61 

0. 75 

0.70 

to predicted values. 

to 

to 

to 

to 

to 

to 


0.88 

1.03 

1. 10 

1. 21 

0.83 

0.88 


(2) Uniform wall temperature; round tubes: Average 
lieat-transfer coefficients in round tubes with uniform wall 
temperatures or with wall conditions somewhere between 
uniform wall temperature and uniform heat input were 
measured by Gilliland, Musser, and Page; and by Doody and 
Younger. In figure 46 are curves representing mean lines 
through the data of Gilliland, Musser, and Page; the data 
of Doody and Younger are represented by a cross-hatched 
area because of scatter. Also shown in figure 46 are the 
relations for average Nusselt number (i/D— 100) calculated 
from equations (1) and (2) and figure 44. The following 
data are not considered: 



Figure 46. — Comparison of measured and predicted average Nusselt 
numbers in round tubes with uniform wall temperature or with a 
wall condition somewhere between uniform wall temperature and 
uniform heat input. 

(a) Data below a Peclet number of 200: These data are 
in the transition flow region. 

(b) Data of Doody and Younger: The scatter is large. 
The Nusselt number of the remaining data can be com- 
pared with the predicted values for an f/D of 100 as follows 
(values of Nusselt number halfway between the values of 
the two theoretical curves in fig. 46 will be used for com- 
parison): 


Peclet number, Pe 

500 

1000 j 

Range of ratio of measured values to predicted values. 

0. 63 

1 

0.61 1 

to 

to 1 


0. 69 

0. 72 j 


(3) Annuli: Average heat-transfer coefficients in annuli 
or between flat plates with constant heat input to the wall 
were measured by Lyon; Werner, King, and Tidball; 
Sineath; Trefethen; and Lubarsky. Figure 47 shows curves 
representing mean lines through the data of these various 
investigators; also shown are the relations for average 
Nusselt number (l/D= 100) calculated from equations (4) 
and (5) and figure 44. The Nusselt numbers of the data 
compare with the predicted Nusselt numbers (average of 
the Nusselt numbers of the two theoretical curves of fig. 
47) as follows (values below a Peclet number of 200 are not 
considered because they fall in the transition flow region): 


Peclet number, Pe 

2C0 

500 

1000 

Range of ratio of measured values to predicted values. 

0.69 

0. 39 

0.39 


to 

to 

to 


1.21 

1.37 

0. 67 
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Figure 47. — Comparison of measured and predicted average Nusselt 
numbers in annuli and between flat plates with uniform heat input 
to wall. 


Temperature distribution. — The only experimental data 
on temperature distribution are those of Isakoff and Drew. 
Plots of the temperature distributions measured by them are 
shown in figure 32; Martinelli’s predicted temperature distri- 
butions are shown for comparison. It is possible to use 
Isakoff and Drew’s temperature and velocity profiles (the 
measured velocity profiles check quite well with the pre- 
dicted velocity profiles) to calculate the values of the ratio 
(t w — t m )/(t w — t c ) shown in figure 48. Martinelli’s predictions 
(fig. 5) for the ratio (t m —t m )/(t w —t c ) are also shown. The 
measured values are smaller than the predicted values. The 
predicted values of Martinelli for (t m —t m )/(t w —t c ) were used 
to calculate the fluid bulk temperature in those cases in 
which fluid centerline temperature was measured (Elser; 
Doody and Younger; Werner, King, and Tidball (test 
section A)). If the values of (t m —t m )/(t m —t c ) are actually 
lower than predicted by Martinelli, the Nusselt numbers of 
these cases would increase. 



Figure 48.— Comparison of measured and predicted values of ratio 
C tv,—t m )l(tw—t c ) in round tubes with uniform heat input to wall. 


1.4 


A Lubarsky, ref. 16 
B Stromquist, ref. 21 

C Johnson Hartnett, and Clabaugh (lead-bismuth 
eutectic}; fully developed data, refs. [7 and 18 
D Isakoff and Drew, fully developed data (inside 
wall temperatures calculated from temper- 
ature profiles), refs. 19 and 20 
E Seban, ref. II 
F Lyon, ref. 6 

G Johnson, Clabaugh and Hartnett (mercury); 

fully developed data, ref. 24 
H Trefethen, refs. 12 and 13 
J English and Barrett refs. 9 and 10 
K Werner, King, and Tidball (test section B), ref. 7 
M Werner. King, and Tidball (test section A) 
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Figure 49. — Variation of ratio of measured Nusselt number to pre- 
dicted Nusselt number with Peclet number. 


Final comparison of heat-transfer data. — The experi- 
mental data of all the investigators for fully developed and 
average heat-transfer coefficients have been compared with 
the appropriate prediction, and the results are shown in 
table I. 

The variation of the ratio of measured to predicted Nusselt 
number with Peclet number is shown in figure 49 for some of 
the data of table I. The results which are not shown in 
figure 49 were not included for the following reasons: 

(1) There is large scatter of data. 

(2) Obvious uncertainties exist as to the accuracy of the 
data. 

(3) The measurements are of average heat-transfer co- 
efficients which were made concurrently with the measure- 
ments of fully developed coefficients; the fully developed 
coefficients are shown in figure 49. 

(4) The measurements are of annulus heat-transfer co- 
efficients which were made concurrently with the measure- 
ments of z-ound-tube coefficients; the round-tube coefficients 
are shown in figure 49. 

On the basis of the results shown in table I and figure 49, it 
can be seen that most of the measured values of fully de- 
veloped and avei'age Nusselt numbei’s for turbulent flow (as 
given by eqs. (1), (2), (4), and (5) and fig. 43) fall between 60 
to 80 percent of their predicted values. 

SUGGESTED EXPERIMENTAL WORK 

It is suggested that the type of experiment most likely to 
reduce the uncertainties with respect to liquid-metal heat 
transfer would be one in which velocity and tempei’ature 
pi’ofiles were measured in the fluid, somewhat like the ex- 
pei'iment of Isakoff and Drew. The experiment of Isakoff 
and Drew could be improved by the use of a thick, high- 
conductivity metallic coating ai’ound the test section similar 
to the one used by Seban or Johnson, Hartnett, and Clabaugh; 
this would probably eliminate the uncertainties in the 
measurements of outside wall temperature. 

The experimental data are insufficient to lead to any con- 
clusion concerning liquid-metal heat ti’ansfer in the laminar 
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and transition flow regions. Such data are greatly needed, 
because the small amount of data in these flow regions 
disagrees considerably with theoretical predictions. 

SUMMARY OF RESULTS 

The review of the experimental investigations of liquid- 
metal heat transfer may be summarized as follows: 

1. The experimental data of the various investigators were 
reevaluated using assumptions and methods as consistent as 
possible, and the results were compared with each other and 
with theoretical values. 

2. The reevaluated experimental data for fully developed 
Nusselt number in the turbulent-flow region were found still 
to have considerable spread, and most of the data are lower 
than predicted theoretically. 

3. An equation based on empirical grounds, which best 
represents most of the fully developed heat-transfer data, is 

Nu= 0.625 Pe 0A 

where Nu and Pe represent Nusselt number and Peclet 
number, respectively. 

4. The theoretical predictions of heat transfer in the 
entrance region were found to give lower values, in most 
cases, than those found in the experimental work. 

5. Integrating the theoretical and experimental results for 
the ratio Nu x /Nu f gave predictions for the value of the ratio 
Nuav/Nu f over a range of Peclet number and length-diameter 
ratio. 

6. The small amount of data on temperature distribution 
disagreed with the theoretical predictions, the discrepancy 
increasing with decreasing Reynolds number. 

7. The experimental evidence is insufficient to lead to any 
conclusion about liquicl-metal heat transfer in the laminar 
and transition flow regions. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 4, 1954 
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TABLE I.— COMPARISON OF HEAT-TRANSFER DATA 


Investigation 

Ref. 

Type of heat-transfer 
coefficient measured 

Theoretical eq. used for 
comparison 

Ratio of measured Nusselt 
number to predicted Nus- 
selt number for Peclet 
number of — 




200 

500 

1000 

2000 

Styrikovich and Semen- 
ovker. 

1 

Round tube, fully devel- 
oped, uniform heat input. 

Nu f = 7.0 + 0.025 Pe t °- 8 

-- 

-- 

0. 80 

0. 77 

Gilliland, Musscr, and 
Page (heating data). 

2,3 

Round tube, over-all av., 
uniform wall tempera- 
ture. 

5. 0 + 0. 025 Pe/- 8 , 
“corrected for l/D = 4 4. 

-- 

0. 69 

0. 72 

-- 

Gilliland, Musser, and 
Page (cooling data). 

2, 3 

Round tube, over-all av., 
between uniform heat 
input and uniform wall 
temperature. 

Nu,= 6.0 + 0.025 Per 0 *, (av. 
of eqs. (1) and (2)), cor- 
rected for l/D— 160. 


0. 63 

j 

0. 61 


Elser__ 

4 

Round tube, fully devel- 
oped, uniform heat input. 

A t h/= 7.0 + 0.025 Pe/- 8 

-- 

"" 

-- 

0. 17 
. 40 

Bailey, Cope, and Wat- 
son. 

5 

Round tube, fully devel- 
oped, between uniform 
heat input and uniform 
wall temperature 

A T »/=6.0 + 0.025 Pc/- 8 , (av. 
of eqs. (1) and (2)). 

0. 36 

0. 4(> 

0. 39 

0. 35 

Lyon (tube data) _ _ _ 

6 

Round tube, over-all av., 
uniform heat input. 

Au/=7.0 + 0.025 Pe/- 8 , cor- 
rected for l/D — 110. 

0. 63 

0. 72 

0. 80 

-- 

1 Lyon (annulus data) 

i 

! 

6 

Annulus, over-all av., uni- 
form heat input. 

A T M/, n7l = 5.8 + 0.020 
Pc/ „ 1J , corrected for 
I/D = 225. 

0. 74 

0. 90 

-- 

-- 

Untermeyer (without 
magnesium additions). 

--- 

Round tube, fully devel- 
oped, uniform heat input. 

Nu f = 7.0 + 0.025 Pe/ 8 


0. 16 

0. 23 

0. 48 

Untermeyer (with mag- 
| nesium additions). 

--- 

Round tube, fully devel- 
oped, uniform heat input. 

Nu f = 7.0 + 0.025 Pe/- 8 


0. 92 

0. 80 

-- 

Werner and King (heat 
; exchanger A, tube 

[ data) . 

— 

Round tube, over-all av., 
uniform heat input. 

A0(/=7.0 + 0.025 Pe/- 8 , cor- 
rected for l/D — 49. 

0. 69 

0. 76 

0. 79 

-- 

! Werner and King (heat 
exchanger A, annulus 
: data). 

i 


Annulus, over-all av., uni- 
form heat input. 

A T «/. a „ = 0.75(D o /D,)»- 3 X 
(7.0 + 0.025 Pe/, o „ 0 - 8 ), 
DJDi— 1.83, corrected 
for l/D = 55. 

0. 90 




Werner, King, and Tid- 
ball (heat exchanger 
B, tube data). 

7 

Round tube, over-all av., 
uniform heat input. 

Nuf= 7.0 + 0.025 Pe/- 8 , cor- 
rected for l/D = 49. 


0. 97 

1. 04 

1. 14 

Werner, King, and Tid- 
' ball (heat exchanger 

B, annulus data). 

7 

Annulus, over-all av., uni- 
form heat input. 

Niif , an = 0.7 5 (D 0 / Di)°- 3 X 
(7.0 + 0.025 Pe/. a „°- 8 ), 
D 0 /Di=\ .83, corrected 
for l/D = 55. 

1. 15 

1. 29 


1 

1 

Sineath - - - ’ 

8 

Rectangular ducts, over-all 
av., uniform heat input. 

Wu/= 5.8 + 0.020 Pe/,„„ 0 - 8 , 
corrected for l/D = 50. 

-- 

0. 41 

0. 40 

-- 

English and Barrett 

9, 10 

Round tube, fully devel- 
oped, uniform heat input. 

Nu/= 7.0 + 0.025 Pe/- 8 

0. 61 

0. 74 

-- 

-- i 

| 

Seban . — 1 

1 

" 1 

Round tube, fully devel- 
oped, uniform heat input. 

A r «/= 7.0 + 0.025 Pe/- 8 

i 

-- i 

-- 

0. 67 

0. 58 j 

Trefethen (fully devel- 
oped tube data). 

12, 13 

Round tube, fully devel- 
oped, uniform heat input. 

Nu f = 7.0 + 0.025 Pe/- 8 

0. 68 

0. 77 

0. 78 

0. 76 

Trefethen (over-all aver- 
age tube data). 

12, 13 

Round tube, over-all av., 
uniform heat input. 

Nut= 7.0 + 0.025 Pe/- 8 , cor- 
rected for l/D = 65. 

0. 74 

1 

0. 83 

0. 86 

0. 87 1 

\ 


» “Corrected for l/D” means that the fully developed Nusselt number found from the eq. was multiplied by the ratio 
Nu a , /Nuf from fig. 43. 
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TABLE I. COMPARISON OF HEAT-TRANSFER DATA— Concluded 


Investigation 

Ref. 

Type of heat-transfer 
coefficient measured 

Theoretical eq. used for 
comparison 

Ratio of measured Nusselt number to 
predicted Nusselt number for Peclet 
number of — - 

200 

500 

1000 

2000 

5000 

9000 

Trefethen (over-all aver- 
age annulus data) . 

12, 13 

Annulus, over-all av., uni- 
form heat input. 

Average of Nu;= 5.8+0.020 
Pe/- 8 and Nu r =0.75 {DJ 
Di) 03 (7.0+0.025 Pe/,an°- s ), 

2, corrected for l/D= 

i 

200. 

0. 81 

0. 87 





Doody and Younger (data 
with no sodium addi- 
tions) . 

14, 15 

Round tube, over-all av., 
between uniform heat 
input and uniform wall 
temperature. 

Nu f = 6.0 + 0.025 Pe/- 8 , (av- 
erage of eqs. (1) and (2)), 
corrected for l/D = 1 1 4. 

0. 22 
. 51 

0. 33 
. 53 

— 


— 

— 

Doody and Younger (data 
with sodium additions). 

14, 15 

Round tube, over-all av., 
between uniform heat 
input and uniform wall 
temperature. 

Aht/= 6.0 + 0.025 Pe/- 8 , (av- 
erage of eqs. (1) and (2)), 
corrected for l/D = 1 1 4. 

0. 50 
. 92 






Lubarsky (tube data) 

16 

Round tube, over-all av., 
uniform heat input. 

Nuf= 7.0 + 0.025 Pc/- 8 , cor- 
rected for l/D— 100. 


-- 

0. 54 

0. 61 

-- 

-- 

Lubarsky (annulus data) 

16 

Annulus, overall-av., uni- 
form heat input. 

Af|(/. a „ = 5.8 + 0.020Pe/. o „ l >- 8 , 
corrected for l/D = 320. 

-- 

0. 59 

0. 72 

-- 


-- 

Johnson, Hartnett, and 
Clabaugh (lead-bis- 
muth eutectic, fully 
developed data). 

17, 18 

Round tube, fully developed, 
uniform heat input. 

Ab</=7.0 + 0.025 Pc/- 8 



0. 76 

0. 70 



Johnson, Hartnett, and 
Clabaugh (lead-bis- 
muth eutectic, over-all 
average data). 

17, 18 

Round tube, over-all av., 
uniform heat input. 

Nu f = 7.0 + 0.025 Pe/- 8 , cor- 
rected for l/D= 74. 

— 


0. 77 

0. 72 



Isakoff and Drew (fully 
developed data, inside 
wall temperature calcu- 
lated from fluid temper- 
ature profile). 

19, 20 

Round tube, fully developed, 
uniform heat input. 

Nu,= 7.0-0.025 Pc/- 8 


__ j 0. 97 
1 

1. 05 

1. 14 

1. 25 

Isakoff and Drew (fully 
developed data, inside 
wall temperature calcu- 
lated from outside wall 
temperature) . 

19, 20 

Round tube, fully devel- 
oped, uniform heat input. 

Nu f = 7.0 + 0.025 Pe/- 8 



0. 95 

0. 94 

0. 90 

0. 91 

Isakoff and Drew (over- 
all average data, inside 
wall temperature calcu- 
lated from outside wall 
temperature). 

19, 20 

Round tube, over-all av., 
uniform heat input. 

Nu/= 7.0+0.025 Pe/- 8 , cor- 
rected for l/D= 138. 



0. 93 

0. 84 

0. 84 

0. 86 

Stromquist _ 

21 

Round tube, fully devel- 
oped, uniform heat input. 

Ah</= 7.0 + 0.025 Pe/- 8 

0. 54 

0. 55 

0. 57 

0. 58 

0. 54 

0. 51 

MacDonald and Quitten- 
ton. 

22, 23 

Round tube, fully devel- 
oped, uniform heat input. 

ATm/= 7.0 + 0.025 Pe/- 8 

1. 0 
. 26 

-- 

-- 

-- 

-- 

-- 

Johnson, Clabaugh, and 
Hartnett (mercury, ful- 
ly developed data). 

24 

Round tube, fully devel- 
oped, uniform heat input. 

IVw/= 7.0+ 0.025 Pe/- 8 

0. 68 

0. 70 

0. 71 

0. 70 

0. 65 

0. 60 

Johnson, Clabaugh, and 
Hartnett (mercury, 
over-all average data) . 

24 

Round tube, over-all av., 
uniform heat input. 

Nuf= 7.0 + 0.025 Pc/- 8 , cor- 
rected for l/D = 74. 

0. 70 

0. 73 

0. 75 

0. 76 

0. 76 

0. 69 

Johnson, Hartnett and 
Clabaugh (laminar and 
transition flow). 

25 

Round tube, fully devel- 
oped, uniform heat input. 

N«/= 7.0 + 0.025 Pe/- 8 

0. 68 
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